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In this talk we present the recent calculation of the NLO QCD corrections to the production of 
four b-quarks induced by initial state quarks at the CERN LHC. We describe the details of the 
calculation using the GOLEM approach for the virtual corrections and MadEvent/MadDipoIe for 
the real emission part and present some results. 
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1. Introduction 



The study of the Higgs sector of the Standard Model (SM) and the search for New Physics at 
the TeV scale at the LHC experiments at CERN will ultimately depend on a detailed understanding 
of signal and background processes involving many particles [§0|,i. 

When considering hadronic collisions the results obtained using only leading order (LO) in 
perturbation theory retain a strong dependence on the unphysical factorisation and renormalisation 
scales which limits the predictivity of the theoretical computations. 

Whereas the inclusion of next-to-leading order (NLO) QCD effects is relatively straightfor- 
ward for simple kinematic situations (2 — )• 2 and 2 — )• 3 processes), processes with four parton final 
states are at the forefront of the present computing capabilities and the number of available results 
is very limited. However a lot of progress has been made during the last couple of years. Recent 
achievements are the evaluation of the order corrections to pp — )• bbtt ^ 0], pp — )• Wjjj 
]g,^,^^,pp^ ttjj [|l|] and pp ^ Z/yjjj 



In this contribution we present the results obtained in [14| for the process pp ^ bbbb which is 
a relevant background for Higgs searches in the minimal supersymmetric extension of the Standard 
Model (MSSM). The MSSM contains two Higgs doublets containing five scalar Higgs bosons. As 
light Higgs bosons prefer to decay into Z^-quarks in large parts of parameter space, experimental 
studies are so far strongly affected by the uncertainty of the Standard Model backgrounds, espe- 
cially pp ^ bbbb, see for example Ref. [|l^]. Although mainly motivated by supersymmetry, the 
ioux-b final state also allows the study of other interesting scenarios for Beyond Standard Model 
(BSM) physics, such as Hidden Valley models, where decays of hadrons of an additional con- 



fining gauge group can produce high multiplicities of bb pairs [g[ [16|]. This is the reason why 
the pp ^ bbbb process was added to the experimenters wishlist of relevant next-to-leading order 
computations 



2. Calculation 



At leading order two partonic subprocesses contribute to the computation of the pp ^ bbbb 
cross-section: the quark initiated subprocess (qq — )• bbbb) and the gluon initiated one (gg — bbbb). 
In this contribution we presents results obtained for the quark induced subprocess Throughout the 
whole calculation we treat the Z^-quarks as massless, an approximation which is well motivated by 
the LHC kinematics and the applied cuts, and neglect the effect of a heavy top quark altogether. 
Furthermore we neglect the ft-quark contribution in the parton distribution functions, due to the 
smallness of the b PDF. 



2.1 Virtual corrections 

The calculation of the virtual corrections is based on the GOLEM approach |[T^,[T8, 1^, 21] for 
one loop amplitudes. It uses a Feynman diagrammatic representation generated by QGRAF [p3|]. 
FORM [24] is used to map the helicity amplitudes to a tensor form factor representation as defined 
in [^, 22, 25]. The latter is exported to Fortran95 code such that it can be linked to the form 
factor library golem9 5 [20]. For the code generation we have developed a dedicated optimisation 
tool which can also be used for other purposes [E3]. To validate the results a second independent 
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calculation has been performed. There we did a symbolic tensor reduction down to genuine scalar 



integrals where we used FeynArts and FeynCalc [ ]27[ ] to generate the amplitudes. The tensor 



integrals are reduced symbolically to scalar integrals also using the formalism described in []17[]. 



More recently we have independently calculated the virtual connections using the SAMURAI 



package [p8[]. SAMURAI is a package for the reduction of one-loop amplitudes at the integrand 
level based on the OPP approach [ ^9| , Combined with a modified version of our matrix ele- 
ment generator the code for the numeric evaluation of the matrix element was generated in a fully 
automated way. 

Using this setup we were able to reduce the CPU time needed for the evaluation of the virtual 
corrections of a single phase space point to around 0.4s, which is around a factor of 10 faster 
compared to the previous approach. 

This improvement in the code speed is a fundamental step towards the computation of the 
gluon induced subprocess which, due to the larger number of Feynman diagrams and their increased 
complexity, is computationally more intensive than the quark initiated one. 

2.2 Real emission 

The tree level and real emission matrix elements have been generated using MadGraph [E^], 



for the subtraction terms we used MadDipole [|29|, |30p which is based on the dipole formalism 



p^ ]. For the full proton initiated process there are three subprocesses contributing: qq — )• 
bbbbg, gg — bbbbg and qg — bbbbq. Here we restrict ourselves to the first one. The subtraction 
terms and the real emission matrix element have been implemented into the MadEvent framework 



[33] 



A second implementation of the subtraction terms in Whizard [ )34| , [35[ ] has been used to ensure 
the correctness of our code. 



3. Results 

In order to provide a check of the numerical stability of our implementation we looked at the 
cancellation of the double and single poles between the virtual corrections and dipole subtraction 
terms. The result is presented in Fig. [I| and shows that the cancellation happens within double 
precision accuracy for both single and double poles, thus confirming that possible numerical insta- 
bilities are well under control. 

The numerical results we present in this section were obtained for = 14 TeV center of mass 
energy. In a first step we cluster the final state particles to Z^-jets using the ^7 -algorithm [|^]. The 
resulting jets have to pass the following cuts: 

pribj) > 30 GeV 
\v{bj)\ < 2.5 

AR{bi,bj) > 0.8 . (3.1) 

One of the main goals of calculating higher orders corrections is to reduce the dependence of the 
result on renormalisation and factorisation scale. In Fig. ^ we compare the dependence on the 
renormahsation scale for the total cross section and for the invariant mass distribution of the two 
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Figure 1: Cancellation of the poles from the virtual corrections and the subtraction terms normalised to the 
Born. To produce this plot we integrated the process qq — >■ bbbb using 10^ phase space points and the cuts 
described in the text. 
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Figure 2: Dependence on the renormalisation scale for the total cross section (left plot) and for the distri- 
bution of the invariant mass of the two fe-jets with the higest pr (right plot) when varying the scale around 
jj.r=x- lio—x - Pribj) with .jc e [ j , 8] (left plot) and [5 , 2] . (right plot). The dashed line in the right 
plot denotes the leading order prediction for x=\/2. 



Z7-jets with the highest pj at LO and NLO. As expected, a clear reduction on the scale dependence 
can be observed. The central scale is chosen to be /Xq = y^Ly Pri^j) ^^id scale is then varied 
with Hr =x- Ho. As we have not included all possible initial states we kept the factorisation scale 
fixed to be /if = 100 GeV. 
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